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The Ising model in a random field and with power-law decaying ferromagnetic bonds is studied at 
zero temperature. Comparing the scaling of the energy contributions of the ferromagnetic domain 
wall flip and of the random field a la Imry-Ma we obtain a threshold value for the power p of the 
long-range interaction, beyond which no critical behavior occurs. The critical threshold value is 
p c = 3/2, at a difference with the zero field model in which p c = 2. This prediction is confirmed 
by numerical computation of the ground states below, at, and above this threshold value. Some 
possible implications for the critical behavior of spin-glasses in a field are conjectured. 
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Introduction. According to the well known Imry- 
Ma argument [l]-[3| the Random Field Ising Model 
(RFIM) with nearest-neighbor interaction does not dis- 
play any spontaneous magnetization in D < 2. Sponta- 
neous magnetization is, instead, present in D = 3 where 
a rigorous result has shown the occurrence of a finite 
dimensional phase transition 0, [E| and numerical sim- 
ulations [6l-llCf confirm this result. Further analysis by 
rigorous approach [ll| . p erturbation theory [l2[ , and RG 
transformations [13, LLJ] have shown that D = 2 is, ac- 
tually, the lower critical dimension, and no evidence has 
been provided for the existence of any transition in di- 
mension two, both for T > and at T = 0. In the latter 
case the relevant variable is the strength of the random 
magnetic held, i.e., the square root of its variance, nor- 
malized to the ferromagnetic coupling. Renormalization 
group arguments show that the finite temperature tran- 
sition is dominated by the zero temperature fixed point. 

In the present work we present the investigation of the 
zero temperature critical behavior in a one-dimensional 
RFIM with long-range (LR) power-law decaying inter- 
action. Our aim is twofold: (a) to characterize the 
threshold value of the power above which the system 
docs not undergo any phase transition; (b) to gain in- 
sight about the correspondence between LR models with 
a certain power of the interaction decay and short-range 
(SR) models in a given dimension D in presence of a 
field. Our main result is that the critical threshold value 
for the power corresponding to a lower critical dimension 
is p c = 1.5 in the ID RFIM, rather than p c = 2 as in the 
ID ferromagnetic model in absence of a field: the latter 
being the well known Kondo problem. [I5l4l7| . This has 



direct consequences on the determination of the lower 
critical dimension in presence of a field by means of the 
analogy between long-range (LR) and short-range (SR) 
systems. 

SR-^LR connection with no field. We recall that a 
quantitative relationship can be established between the 
power-law p of the LR interaction decay in a ID lattice 
and the dimension D of a SR system displaying the same 



critical behavior. The requirement that the renormalized 
coupling constant has the same scaling dimension leads 
to 



(1) 



Below the Upper Critical Dimension (UCD), though, i.e., 
for p > p m f [40|| , such relationship is not exact anymore. 
Moreover, it grossly fails at the Lower Critical Dimen- 
sion (LCD), D = 1 for the purely ferromagnetic model, 
predicting a p c = 3 > 2 in a ID LR chain. We can im- 
prove Eq. ([1]) by looking at the behavior of the renormal- 
ized space correlation function at criticality in SR model: 
C(r) ~ r -D+2-r, sr (D) _ R cqu i r i n g that at the LCD the 
correlation function does not display any power-law crit- 
ical decay, i.e., D = 2 — r] sr (D) and imposing the correct 
p c = 2, Eq. (fT]) is modified as 
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By construction it is exact at the LCD. The same relation 
holds for Heisenberg ferromagnets, where at the LCD 
(D = 2), r] SI (D) = 0. Eq. © has been first obtained, in 
the framework of spin-glasses, by comparing the singular 
part of the free energy per spin in a LR system of N = L 
spins and in a D-dimensional SR system with the same 
number of spins, N = LP . The magnetic scaling expo- 
nents turn out to follow the relationship y\ = y^'(D)/D, 
being 2yh = D + 2 — 77 [3, Ell ■ Since in LR models, both 
with and without quenched disorder, the two point vertex 
function is not renormalized and rji r = 3 — p [20T - 22 1 also 
in the infrared divergence regime, Eq. ([2]) is recovered 

3 El mull- 

Eq. © states that the critical behavior of the two 
models, i.e., the D-dim. SR and the ID LR mod- 
els, should be similar for all (p, D) couples between 
(p mU UCD) and ( Pc ,LCD) gl|. For p < p ml = 3/2 
the system is in the mean-field regime. In the ID Ising 
model without field this corresponds to D > Z?ucd = 4. 
As D < Ducd infrared divergences occur in the vertex 
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functions and a non-zero anomalous exponent. In D = 3 
a good numerical estimate is r\ = 0.031(5) 25J, corre- 
sponding to p = 1.656(2). In D = 2, Onsager solution 
yields ?y sr = 1/4 and the system is "critically equivalent" 
to p = 15/8 LR model. 

By direct inspection, it is known that no transition is 
present at p > p c = 2. Exactly at p = 2, though, a phase 
transition docs occur. This is the Kondo transition in ID 



magnetic chains [17[. On the contrary, the SR ID Ising 
chain does not display any critical point. This is, actu- 
ally, not unusual and it is due to a direct long interaction 
of iterfaces in LR models. The critical behavior of the LR 
model at p c and of the SR model exactly at the LCD is 
often different: in some instances no transition is present 
in the SR model, while a transition may be present in 
the corresponding LR model. We anticipate that this is 
the case for the RFIM as well: no transition at the SR 
LCD, but a T = fixed point with logarithmic scaling in 
the LR model at p c . 

We stress once again that the same Eq. holds for 
systems with quenched bond disorder, the so-called spin- 
glasses, in which a rigorous result confirms p c = 2 [261 ] . 
Only the mean-field threshold value of p is modified, be- 
cause the relevant interaction term at criticality, and, 
thus, the upper critical dimension (UCD), is different: 
PZ =4/3 E" ' 
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SR^-LR connection in a field. As an external field 
is switched on a new critical fixed point arises that is dif- 
ferent from the zero-field fixed point. This is true both 
for systems with and without quenched bond disorder. 
Lower and upper critical dimensions appear not to de- 
crease in all known cases. In particular, the critical di- 
mensions of the RFIM increase to become Djjcd = 6 
and Dlcd = 2. The extension of Eq. [2] to the random 
magnetic case requires some care. Different definitions 
of the exponent ?7 sr are, indeed, possible since connected 
and disconnected correlation functions decay differently 
and hyper-scaling does not hold [fHH, 27-29j. Here, we 
define an exponent fj sv by the condition that the Fourier 
transform of spin-spin disconnected correlation behaves 
in momentum space as fc _4+ '' Br , or equivalently in posi- 
tion space C^ r ISC (r-) ~ r -- D + 4 ~'fer(.D) ; w h ere the Schwartz- 
Soffer inequality holds: fj sr < 2?y sr J28[. The difference 
between 2ij and fj decreases with the dimension 0, [l(| , 
eventually tending to zero at the LCD. 

We, now, present our study of the ID LR RFIM. First, 
using an Imry-Ma-like argument we predict p c = 1.5. 
Further, we analyze the critical behavior of the model at 
p ~ p c by means of numerical computations of the ground 
states properties at T = as function of the strength of 
the ferromagnetic interaction J. 

The long-range RFIM and the Imry-Ma argument. 
The Hamiltonian of the LR ID RFIM is 



where Jij = J\i — j\~ p and hi is a random field with a 
bimodal distribution of zero average and variance h . 

In a ordinary fcrromagnet the cost to flip a domain 
of spins of length L grows like L 2 ~ p . As the random 
field is switched on this will compete with the energy of 
the orientation along the field going like L 1 ! 2 . Accord- 
ing to the argument developed by Imry and Ma for SR 
D-dimensional systems [l|, as p > 1.5 there will always 
be a size large enough for the field to destroy any ferro- 
magnetic domain and no long-range order can be estab- 
lished. The exponent value p = 1.5 should, therefore, be 
the analogue of the LCD in nearest-neighbor interacting 
L>-dimcnsional RFIM, i.e. D = 2 [3-11, i, [Io| . 

Levy lattice. In order to validate this analytic pre- 
diction we performed numerical estimates of the ground 
state properties at T = for the ID RFIM on a Levy 
lattice, that is a finite connectivity random graph equiv- 
alent to a fully connected LR model [3(|. In this dilute 
graph two sites i and j are connected (i.e., Jij ^ 0) with 
a probability 



P(Jij = J) 



(4) 



where the sum runs over all possible distances realizable 
on the ID chain of length L and such that the total num- 
ber of bonds is independent from p and equal to zL, 
where z is the average spin connectivity. For p large 
enough one has a nearest-neighbor chain, whereas for 
p = the distribution of the connectivities is Poissonian 
and the system corresponds to an Erdos-Rcnyi graph. 

Numerical results. Using the Minimum Cut algo- 
rithm 0, [3l], [32| of the Lemon Graph Library [33| we have 
computed thermodynamic observables on T = ground 
states of Levy graphs of different length, averaging over 
different realizations of random fields. The computation 
has been performed varying the ferromagnetic coupling 
magnitude J. 

First, we present the numerical results at p = p c = 1.5 
where we used sizes ranging from L = 250 to L = 256000. 
The number of samples of disorder are 10000 for L < 
64000, 5000 at L = 128000 and 2000 at L = 256000. For 
each sample we compute the ground state for 41 values 
of the J coupling in the interval [0.2 : 0.4]. The random 
field mean square displacement is kept constant, h = 1. 

To understand whether a critical behavior is there we 
study the finite size behavior of the Binder cumulant 
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U 



(3) 



If, in the thermodynamic limit, a phase transition occurs 
at a given critical field h c , the Binder cumulant will be 
one (long-range order) for h < h c and zero for h > h c . 
As L increases, we observe that the various Binder curves 
tend to a limiting curve, cf. Fig. [TJ with a behavior that 
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FIG. 1: Finite size Binder cumulants versus the strength of 
the random field in J units at p — 1.5. The critical value 
estimate by FSS analysis is (h/J) c = 2.31(5). 
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FIG. 3: Rescaled Binder cumulant versus ft/ J at p — 1.5. 
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FIG. 2: Multilinear interpolation for the J/h value at which 
finite size Binder cumulants significantly change value (g = 
0.3, 0.4, 0.5, 0.6 and 0.7) with J L (]nL\g) = J c + b g /ln L. 



we will show compatible with a scaling logarithmic de- 
cay. To estimate the critical point we look at the values 
of h/J that, for different sizes, yield the same g value. 
Specifically, in Fig. [2] we show the behavior of the lim- 
iting inverse critical field value, (J/h) c , as computed at 
every size for different fixed values of the Binder cumulant 
(g = 0.3, 0.4, 0.5, 0.6, 0.7) versus (log L)~ l . In the L — > oc 
limit all curves are compatible with a multiple linear fit 
in 1/lni yielding the estimate (J/h) c = 0.433(9), or 
(h/J) c = 2.31(5). 

In Fig. [3]we plot the curves in the rescaled variable and 
observe a very good overlap in the critical region. To fur- 
ther characterize the transition we look at the behavior 
of magnetization momenta at the critical point. In Fig. 
E] we, thus, present the behavior of the squared magne- 
tization around the estimated critical value compatible 
with a logarithmic finite size scaling (FSS). 
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FIG. 4: Rescaled m2 = (s) 2 curves versus h/J at p — 1.5. 



As a comparison, we also present the behavior of the 
Binder cumulant for values of the power p slightly be- 
low and above p c . For p = 1.4 and p = 1.6 we compute 
the ground states of systems of size between L = 250 
and L = 128000 averaging over 10000 disordered fields 
configurations on 51 J values. At p — 1.4, cf. Fig. 
Binder curves cross each other at finite h/J and a 
FSS analysis of the crossing points yields a critical value 
(h/J) c = 3.23(7). Using the scaling property of the 
disconnected correlation function (s) 2 ~ L 3 "^ 11 , with 
fjlr = 3(2 — p) (this formula should hold for p G [1,3/2], 
from the FSS of the crossing points (cf. inset of Fig. 
[5]) we further obtain the estimate (h/J) c = 3.266(2) 
and from the FSS of the derivatives of (s) 2 we estimate 
l/v = 0.316(9). 

For p = 1.6, on the contrary, no crossing is observable 
and the non-zero Binder values continuously run away 
towards larger and larger fields, cf. Fig. El compatibly 
with the claim of absence of transition above p = 1.5. 

Eq. @ is a not too bad approximation for what con- 
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FIG. 5: Finite size Binder cumulants at p = 1.4 for L — 
250, 128000. The critical field estimate is (h/J) c = 
3.23(7). Inset: scale invariant m 2 L ' 2 = (spL*'" 3 vs. h/J, 
77 = 6 - 2p. 

cerns the transition without field. With no field, at 
p = 3/2, Eq. ([2]) would predict a mean-field transi- 
tion (corresponding to D = 4); non- mean- field transi- 
tions would be expected at p = 1.6545 (corresponding 
to D = 3, with ?/ sr (3) = 0.0364(5) [H), and p = 1.875 
(corresponding to D = 2, r; sr (2) = 1/4 [35|); eventu- 
ally, the "LCD" -equivalent exponent value corresponding 
to D = 1 (r? S r(l) = 1) woud be p c = 2. 

Such predictions have been recently numerically inves- 
tigated showing that the critical exponents at p = 1.6546 
and 1.875 do not strictly correspond to, respectively, 2D 
and 3D critical exponents |3(|: if for 3D, nearer to the 
mean-field threshold, numerical estimates are still consis- 
tent with each other, in 2D they appear not compatible 
anymore. A similar trend has been identified in spin- 
glasses where LR systems with values of p equivalent to 
4D and 3D have been analyzed flij ]. 

When the random field is switched on and a new fixed 
point for the RG flow arises the situation changes. The 
mean- field threshold is now p m { = 4/3 (UCD=6). The 
Imry-Ma argument and the simulations presented in this 
work clearly show that the threshold for the critical be- 
havior is p^ = 1.5. 

Discussion and conclusions. From the present 
work we clearly understand that the reference values for 
p are different from those obtained in absence of a field. 
We obtain such evidence by means of an Imry-Ma-like 
argument and a numerical study of the zero temperature 
ground states of the RFIM on a Levy lattice for system 
sizes ranging fro 250 to 256000 spins. Specifically, we find 
that in the ID RFIM with LR interactions no transition 
is present for p > 1.5 and that at p = 1.5 a T = fixed 
point is still present with a logarithmic scaling. 

In presence of a random field we can reformulate the 
"/? — D" relationship Eq. © in terms of the anoma- 



FIG. 6: Finite size Binder cumulants at p = 1.6 for L = 
250, . . . , 64000. 

lous exponent fj sr (D), rather than rj ST (D). That is, we 
consider the most divergent correlation function at criti- 
cality in a SR system in dimension D: the disconnected 
one. At the lower critical dimension (D — 2), where 
D ~ A + fj sr (D) ~ 0, the threshold value of the power p 
has to be equal to the maximum one compatible with the 
existence of a transition: p c = 3/2. This leads to 

,-l-m&>* (6) 

yielding the value of p corresponding to a SR model in 
D dimensions. As Eq. @ in zero field, Eq. is exact, 
at all events, at D =UCD and LCD. Since in the latter 
case fj sr ~ 2rj sr we notice that in this particular case Eq. 
© coincides with Eq. For both Eqs. ^ and © the 
LCD equivalent value of p is the correct one: p c = 1.5. It 
is important to stress that a given value of p corresponds 
to completely different critical behaviors and to different 
dimensions of short-range critically equivalent systems 
if the field is present or absent. As an instance p = 
1.5 is the mean- field threshold in the Ising ferromagnetic 
model, corresponding to UCD D = 4 and it is the critical 
threshold in the RFIM, corresponding to LCD D = 2. 

Does this relationship hold also in presence of random 
bonds, besides random fields ? The Imry-Ma argument 
is specific for the RFIM and cannot be exported to spin- 
glasses because these more complicated systems lack any 
long-range order in the frozen phase. Therefore, a quan- 
titative estimate of the threshold value p^ corresponding 
to the SR LCD is beyond the reach of the analysis pre- 
sented here. Contrarily to the ordered bond model, in 
spin-glasses the UCD does not increase by applying a 
random field (p m f = 4/3). Nor the LCD, that remains 
equal to D = 2.5 according to a computation of interface 
free energy j3?| ■ 

In LR systems in a field, though, we would be surprised 
to observe a spin-glass phase for values of p > pj, 1 = 3/2 
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for which no ferromagnetic transition is present in ab- 
sence of bond disorder. This suggests that much caution 
should be taken in numerical data interpretation in LR 
spin-glass systems in presence of a field. Above all when 
ultrametricity 38( or lack of Almeida- Thouless transition 
}39j are tested at p > 1.5. 
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